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1.  INTRODUCTION 


Irprep  is  a  collection  of  eight  programs.  Three  extract  information  from  the  BRL-CAD  (Ballistic  Research 
Laboratory  Computer  Aided  Design)  three  dimensional  models  and  create  input  files  for  use  in  different 
thermal  models.  The  other  five  programs  are  used  to  display  the  output  from  thermal  models.  The  BRL- 
CAD  models  are  developed  using  the  multi-device  graphics  editor  (mged)  and  use  three  dimensional 
combinatorial  solid  geometry  (CSG).  These  models  are  built  by  combining  different  primitives  (e.g., 
spheres,  cylinders,  and  boxes)  using  boolean  operations.  These  models  are  then  evaluated  using  ray-tracing 
programs  that  fire  rays  at  the  target  and  extract  geometric  data  at  each  intersection  of  the  rays  with  the 
primitives  in  the  model  [1,2]. 

The  eight  different  programs  in  irprep  are  firpass,  seepass,  shapefact,  display,  see,  8ee2,  pictx,  and 
pictsgi.  All  eight  are  interactive  programs.  Firpass  must  be  run  before  seepass,  shapefact,  and 
display;  since  firpass  creates  files  used  by  these  programs.  Display  creates  a  file  used  by  see  and  see2. 

Firpass,  seepass,  and  shapefact  are  the  three  programs  in  irprep  that  are  designed  to  extract  geometric 
data  from  BRL-CAD  models.  Firpass  creates  several  different  types  of  files:  a  diagnostic  file,  an  error  file, 
two  intermediate  files,  and  a  primary  file.  There  are  three  different  types  of  primary  files  that  firpass  may 
create;  the  facet  records  for  PRISM  (Physically  Reasonable  Infrared  Signature  Model)  [3,4],  a  generic  file, 
and  a  geometric  properties  file.  PRISM  is  a  thermal  model  developed  at  Keweenaw  Research  Center  (KRC) 
under  contract  to  the  Tank  Automotive  Command  (TACOM).  Seepass  creates  a  diagnostic  file,  an  error 
file,  an  intermediate  file,  and  a  primary  file.  The  primary  file  is  either  the  conductivity  table  for  PRISM  or 
a  generic  file.  Shapefact  finds  the  shape  factors  between  regions  in  the  engine  and  creates  one  of  two  types 
of  primary  files,  a  regular  or  generic  file  as  well  as  an  error  file.  These  shape  factors  may  be  used  with  the 
engine  compartment  radiation  models  developed  by  Jim  Rapp  at  the  Ballistic  Research  Laboratory  (BRL) 
[5]  or  in  conjunction  with  other  thermal  models. 

The  other  five  programs  in  irprep  are  used  for  post  processing.  Display  takes  a  temperature  output  file 
(either  from  PRISM  or  a  generic  file),  the  three  dimensional  geometry  file,  and  an  intermediate  file  from 
firpass  and  creates  an  output  file  that  can  be  displayed  by  Analyzer2  [6,7],  see,  or  see2.  Display  uses 
ray-tracing  and  the  temperature  output  file  from  a  thermal  model  to  create  temperature  maps.  Analyzer2  is 
a  program  to  display  temperature  information  written  by  Glenn  Durfee  a  Science  and  Engineering 
.Apprentice  Program  (SEAP)  student  with  the  System  Engineering  and  Concepts  Analysis  Division 
(SLC.AD).  Analyzer2  uses  X-VVindows  with  Motif.  See  and  8ce2  both  produce  a  color  display  of  the  target 
where  the  color  represents  the  temperature  of  the  region.  See  displays  the  image  using  X-Windows,  see2 
displays  the  image  using  Silicon  Graphics  (SGI)  graphics.  Pictx  and  pictsgi  combine  the  programs  display 
and  see  or  display  and  see2  respectively. 


2.  MODIFICATIONS  TO  GEOMETRY  FILE 

If  the  geometry  file  being  used  with  firpass  and  seepass  was  not  built  specifically  for  heat  flow  analysis  it 
will  probably  needs  to  be  modified.  Geometric  models  may  be  built  for  any  number  of  purposes: 
vulnerability  evaluation,  radar  signature  analysis,  optical  signatures  determination,  and  infrared  (IR) 
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signature  simulation.  For  each  of  these  uses,  the  models  may  be  built  differently  and  must  be  adapted  for 
an  unintended  use.  Any  regions  that  do  not  have  the  same  thermal  properties  over  the  entire  region  will 
need  to  be  changed.  These  regions  must  be  split  into  smaller  regions  such  that  each  region  has  the  same 
thermal  properties  over  the  entire  region.  Similarly,  any  large  region  may  require  subdivision  into  smaller 
regions  if  large  temperature  variations  can  occur  within  the  region.  Also,  any  region  having  more  than  one 
surface  normal  needs  to  be  broken  into  smaller  regions.  For  example,  if  the  geometric  model  of  a  tank 
turret  is  one  region  or  even  several  large  regions  it  probably  needs  to  be  broken  into  smaller  regions  so  that 
there  is  more  than  one  normal  associated  with  the  turret.  PRISM  needs  one  normal  per  region  and  if  a  large 
region  is  broken  into  smaller  regions  the  normal  will  be  averaged  over  a  smaller  area.  Another  example  is 
toolboxes  made  from  one  solid  or  two  solids  (where  one  is  the  inside  subtracted  from  the  other).  Assuming 
these  boxes  are  placed  with  the  bottom  on  a  surface  they  will  have  a  surface  normal  pointing  straight  up. 
Surface  normals  are  used  in  PRISM  for  solar  loading,  and  appropriate  solar  loading  is  needed  on  the  sides  as 
well  as  the  top  of  the  box.  It  might  be  appropriate  to  model  each  of  these  sides  as  a  separate  region  so  that 
each  side  has  its  own  normal  vector. 

When  firpass  finds  the  centroid  of  a  region  it  will  not  necessarily  be  within  that  region,  for  example,  a 
turret  ring.  A  turret  ring  is  often  made  by  subtracting  two  concentric  cylinders  from  each  other.  Firpass 
will  calculate  the  centroid  as  being  at  the  center  of  the  cylinders  not  taking  into  account  that  the  center  is 
not  part  of  the  region.  If  the  user  does  not  want  the  centroid  to  be  outside  of  the  region,  some  regions  will 
need  to  be  remodeled. 

Any  region  that  is  almost  totally  surrounded  by  exterior  air,  such  as  the  gun  tube  will  need  to  be  examined. 
The  average  normal  vector  to  the  exterior  surface  of  a  gun  tube  should  be  (0,  0,  0)  (where  this  notation 
denotes  vector  components  in  the  x,  y,  and  z  directions).  Firpass  will  do  unpredictable  things  to  a  zero 
normal  vector  when  normalizing  it.  What  happens  is  that  the  normals  do  not  cancel  each  other  out 
completely  and  when  a  vector  is  normalized  it’s  length  must  be  equal  to  one.  One  possibility  for  solving  this 
problem  is  to  split  the  gun  tube  cylinder  into  an  inside  and  outside  cylinder  and  then  splitting  these 
cylinders  lengthwise  into  quarters  or  eighths.  If  there  is  just  one  cylinder  split  into  quarters  there  is  still  the 
problem  that  there  is  exterior  air  on  both  sides  of  the  guntube  (on  the  outer  surface  and  on  the  interior  of 
the  barrel)  and  it  would  be  expected  that  the  normal  would  be  near  (0,  0,  0).  Or  the  user  could  also 
designate  the  air  inside  the  gun  barrel  to  be  interior  air. 

Any  region  not  attached  to  something  will  need  to  be  modified.  Often  for  vulnerability  modeling  toolboxes 
will  not  get  placed  exactly  on  a  fender,  instead  they  will  be  floating  slightly  above  the  fender.  These 
Tolboxes  will  need  to  be  moved  down  or  exterior  air  will  be  detected  below  them.  Any  interior  components 
not  necessary  for  thermal  modeling  may  be  stripped  out,  such  as,  crew,  ammunition,  or  computers.  The 
engine  components  should  be  left  in  the  file. 

Appropriate  air  must  also  be  placed  in  the  model.  When  irprep  finds  the  exterior  surface  area  of  a  region, 
any  surface  bounded  by  exterior  air  or  a  gap  will  be  considered  an  exterior  surface.  For  example,  if  a  box  is 
hollowed  out  and  no  air  is  put  in  this  hollow  space  the  irprep  programs  will  assume  this  is  exterior  air.  In 
the  mged  geometry  file  an  air  region  must  have  a  region  identification  number  (id)  of  zero  and  an  air  code 
greater  than  zero.  The  region  id  indicates  what  type  of  material  the  region  is  composed  of  when  used  with  a 
material  data  base.  The  region  id  and  air  code  are  defined  using  mged  and  may  be  checked  by  listing  the 
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region  (1  region  name)  while  using  the  mged  editor.  They  may  be  changed  by  using  regdef  or  edcodes.  The 
following  air  codes  are  used  in  irprep. 

1  =>  exterior  air 

2  =  >  crew  compartment  air 

5  =  >  engine  compartment  air 

6  =  >  closed  compartment  air 

7  =  >  exhaust  air 

8  =  >  generic  air  1 

9  =>  generic  air  2 

Closed  compartment  air  is  used  to  include  air  in  toolboxes,  this  is  air  that  is  interior  air  but  not  crew 
compartment  air,  engine  compartment  air,  or  exhaust  air.  Generic  airs  one  and  two  have  been  added  for 
future  use  or  if  the  user  is  interested  is  areas  surrounded  by  other  air.  It  has  also  been  suggested  that  the 
material  and  line  of  sight  for  an  air  region  be  set  to  zero. 

One  suggestion  when  building  these  three  dimensional  models  is  to  group  all  groups  except  air  in  one  big 
group  (such  as  all).  All  air  should  be  in  one  group  (such  as  air).  The  model  must  have  air  in  it  for  firpass 
and  secpass  to  run  correctly.  PRISM  does  not  want  air  in  the  vehicle  file.  If  the  groups  and  regions  are 
grouped  this  way  it  is  much  easier  to  delete  unwanted  air  regions  from  the  PRISM  vehicle  file.  This  will  be 
discussed  in  more  detail  later. 

The  geometry  to  be  used  with  shapefact  should  be  the  engine  compartment  consisting  of  engine 
components  and  surrounding  walls  only.  All  interior  air  in  this  part  of  the  description  must  be  modeled  as 
engine  air.  There  must  be  engine  air  in  this  model  for  shapefact  to  work  correctly. 

The  mged  file  that  is  used  as  an  example  throughout  this  user’s  manual  is  test.veh.g.  This  is  a  very  simple 
vehicle  (21  regions  including  three  air  regions).  See  Figure  1  for  an  optical  rendering  of  this  vehicle.  This 
geometry  file  is  not  appropriate  for  thermal  modeling.  It  contains  some  long  plates,  for  example,  the  side  of 
the  hull  is  six  meters  long.  If  this  entire  side  does  not  have  the  same  thermal  properties  over  it’s  entire 
length  it  should  be  broken  into  smaller  regions.  A  wheel  has  also  been  designed  as  being  one  region,  a 
cylinder.  The  normal  vector  of  a  cylinder  is  (0,  0,  0);  as  mentioned  before  firpass  will  do  unpredictable 
things  when  a  region  has  a  normal  vector  of  (0,  0,  0).  The  wheel  should  be  split  into  smaller  pieces.  Regions 
for  this  example  were  not  subdivided  in  order  to  keep  the  number  of  regions  small. 


3  FIRPASS 

3.1.  Firpass  Input  Files 

Firpass  interrogates  an  mged  file  to  obtain  the  facet  records  for  either  a  PRISM-ready  file,  a  generic  file,  or 
a  geometric  properties  rile.  In  order  to  run  firpass  the  user  will  need  a  compiled  copy  of  firpass,  an  mged 
geometry  file,  and  a  material  properties  file.  The  geometry  file  should  have  all  the  necessary  changes  made 
to  it  (see  section  2). 
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Figure  1.  Optical  Rendering  of  Simple  Vehicle 
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The  material  properties  file  is  a  list  of  forty-one  different  materials  that  could  be  contained  in  the  mged  file. 
There  are  seven  columns  in  the  file.  Each  column  is  described  below. 

column  1  -  material  code 

column  2  -  density  (kg/m  ) 

column  3  -  specific  heat  (J/kgC) 

column  4  -  absorptivity 

column  5  -  emissivity 

column  6  -  thermal  conductivHy  (W/mK) 

column  7  -  description  ot  material 

Material  zero  must  be  air  and  entries  for  columns  two  through  six  will  be  zero.  The  last  material  must  be 
number  forty  and  all  materials  between  zero  and  forty  must  exist.  They  may  be  filled  in  with  .zeroes.  A 
sample  material  properties  file  is  located  in  appendix  A.  In  this  material  properties  file  absorptivity  and 
emissivity  have  been  set  to  a  default  value  (1.0).  Since  these  properties  depend  on  the  surface  coating  of  the 
material,  they  should  be  changed  to  appropriate  values  for  specific  models.  For  our  example  absorptivity 
was  set  to  0.60  and  emissivity  was  set  to  0.94. 


3.2.  Firpass  Output  Files 

There  will  be  four  or  five  output  files  generated  when  running  firpass:  a  diagnostic  file  (optional),  an  error 
file,  two  intermediate  files  (to  be  used  with  secpass,  shapefact,  and  display),  and  a  primary  file.  The 
diagnostic  file  is  an  optional  file  used  mainly  for  debugging  purposes.  The  user  may  choose  to  have  this 
information  written  to  a  file,  the  screen,  or  not  at  all.  The  average  user  probably  does  not  need  to  see  this 
information,  it  contains  intermediate  steps  in  getting  the  information  he  is  interested  in  finding. 

The  error  file  will  contain  records  describing  any  region  overlaps  found  in  the  course  of  ray-tracing  the 
geometric  model  and  a  list  of  the  number  of  adjacent  regions  for  each  region.  Provision  is  also  made  for 
listing  to  this  file  any  large  variances  between  surface  areas  and  volumes  of  regions  when  these  quantities  are 
redundantly  calculated  by  firing  rays  in  different  directions.  A  planned  future  revision  of  the  firpass  code 
would  include  the  capability  to  fire  sets  of  rays  in  three  orthogonal  directions,  whereas  in  the  current  version 
only  one  set  of  rays  are  fired  from  a  single  direction. 

The  intermediate  file  created  for  use  with  secpass  contains  the  total  number  of  regions  in  the  model.  For 
each  region  in  the  model  the  file  contains  the  centroid  of  the  region,  material  identification,  and  the  shared 
surface  area  for  each  adjacent  regioii.  The  intermediate  file  created  for  use  with  shapefact  and  display 
contains  the  region  number  and  name.  This  file  allows  the  engine  components  to  have  the  same  region 
numbers  as  the  corresponding  items  found  with  firpass  and  secpass.  This  file  also  makes  it  possible  to 
look  at  a  PRISM  temperature  output  file  and  have  a  region  number  correspond  to  a  region  name. 

The  user  may  create  any  one  of  three  different  primary  files:  facet  records  for  PRISM,  a  generic  file,  or  a 
geometric  properties  file.  The  PRISM  file  will  contain  the  facet  records  needed  by  PRISM.  When  this  file  is 
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combined  with  the  conductivity  table  created  by  secpass  it  becomes  a  PRISM-ready  vehicle  data  file.  The 
information  in  the  PRISM  file  for  each  region  is  listed  below. 

region  number 
region  name 
region  type 

exterior  surface  area  (ra^) 
mass  (kg) 

specific  heat  (J/kgC) 
emissivity 

surface  normal  unit  vector 
convection  coefficient 
region  number  seen  by  back  surface 
region  number  seen  by  front  surface 
shape  factor 
hub  radius 

bearing  friction  constant 

The  facet  records  for  PRISM  are  written  in  the  format  specified  for  PRISM’s  vehicle  data  file.  (See  the 
PRISM  User’s  Manual  [3,4]  for  the  format  that  this  file  is  written.)  The  PRISM  file  may  be  written  so  that 
it  can  be  read  by  PRISM  2.0  or  PRISM  3.0.  There  are  some  very  slight  differences  between  the  files  read  by 
these  two  releases  [3,4]. 

A  number  is  arbitrarily  assigned  by  the  ray  tracer  to  each  region  and  one  is  added  to  it  to  get  the  PRISM 
number.  (PRISM  numbers  must  start  at  one  and  the  numbers  assigned  by  the  ray  tracer  start  at  zero.)  The 
name  is  the  twenty-five  right-most  characters  of  the  region  path  name  [ij.  PRISM  only  accepts  twenty-five 
characters  but  the  region  path  name  may  be  greater  than  twenty-five.  The  surface  normal  unit  vector  is  the 
normalized  average  of  the  individual  normals  computed  at  the  intersections  of  rays  with  the  exterior  free 
surface  area  of  the  region.  The  exterior  surface  area  is  computed  by  adding  over  the  entire  exterior  surface 
area  the  small  area  associated  with  each  ray  that  hits  the  region.  The  volume  is  computed  by  adding  over 
the  entire  region  the  volume  tube  (rectangular  tube)  associated  with  each  ray  that  hits  the  region.  Once  the 
volume  of  a  region  is  known,  its  mass  may  be  found  by  multiplying  its  volume  by  its  density  (from  the 
material  properties  file).  The  specific  heat  and  emissivity  are  dependent  on  the  material  type  and  are 
assigned  from  the  material  properties  file.  The  hub  radius  will  be  set  to  zero  unless  that  region  is  bounded 
by  engine  air  in  which  case  the  area  bounded  by  engine  air  will  be  assigned  to  this  data  field,  (this  field  is 
available  because  the  radius  is  only  used  for  suspension  components).  Region  type  and  convection  coefficient 
are  both  set  to  one;  region  number  seen  by  the  back  and  front  surfaces,  shape  factor,  and  bearing  friction 
constant  are  set  to  zero.  These  variables  are  only  used  for  special  cases  and  it  is  up  to  the  user  to  change 
them  where  necessary.  See  the  PRISM  User’s  Manual  [3,4]  for  a  description  of  when  the  special  conditions 
prevail. 

The  user  may  choose  the  second  option  for  the  primary  file  to  create  a  generic  file  that  has  the  required 
information  for  use  in  other  thermal  models.  The  user  then  can  re-format  the  information  in  the  desired 
format  for  use  in  another  thermal  model.  The  information  contained  in  this  file  for  each  region  is  listed 
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below. 


3  3 

1  region-number  name  volume_{m  )  density _(kg/m  )  thermal-conductivity _(W/mK) 

specific-heat-(  J/kgC )  region-material 

2  region— number  number— of— exterior-convection— surfaces  number_of— interior_convection_surfaces 
area— of-all_exterior-convection_surfaces-(  m  ) 

area— of— all— interior_convection-Surfaces-(m  ) 

3  region— number  number-of_solar-loaded_surfaces 

2 

surface-area_(m  )  surface-norm^’_components  absorptivity 

The  first  digit  of  each  main  line  indicates  the  type  of  information  in  the  succeeding  line.  1  indicates 
geometric  and  material  properties  of  the  region,  all  of  which  is  contained  on  one  line.  A  2  indicates 
convection  heat  transfer  information.  There  may  be  one  to  three  lines  of  this  information.  The  digit  3 
indicates  solar  heat  transfer  and  there  will  be  two  lines  to  contain  this  information. 

The  user  can,  for  a  third  choice,  create  a  geometric  properties  file.  Two  formats  are  possible  for  this  file.  In 
the  first  format,  the  data  are  spaced  so  that  they  are  easily  readable  by  the  user.  In  the  second  format,  the 
data  are  pushed  together  and  there  are  no  columns.  This  format  provides  a  data  file  that  is  easily  read 
using  a  computer  program  for  use  elsewhere.  The  first  format  is  presented  below. 

region-number  name  x-component_of-centroid-(m)  y-component-of_centroid_(m) 

z-component-of_centroid-(m)  volume_(m  )  mass-(kg)  (both  lines  appear  on  same  line) 

region-number  external_surface_area_(m^)  engine_surface_area_(m  )  crew_surface_area_(m*) 

closed_compartment_surface_area-(m  )  exhaust_surface_area_(m  )  generic_l-surface-area_(m  ) 

2 

generic_2_surface_area-(m  )  (all  lines  on  one) 

region-number  material-code  density-(kg/m  )  specific_heat_(J/kgC)  absorptivity  emissivity 
thermal_conductivity_(W/mK)  material-description  (both  lines  on  one) 
region-number  alLadjacent-regions  (separated  by  commas) 

A  maximum  of  twenty  adjacent  regions  may  be  listed. 

All  of  the  data  in  the  second  format  are  in  a  single  string  separated  by  commas  only.  For  each  region,  the 
following  information  is  printed. 

region-number,  name,  x-component_of-centroid_(m),  y-component_of-centroid-(m),  z-component- 

3.2  2 

of_centroid_(m),  volume_(m  ),  mass_(kg),  exterior_surface_area— (m  ),  crew-surface— area-(m  ), 

2  2  2 
engine— surface_area-(m  ),  closed_compartment_surface-area— (m  ),  exhaust-surface_area_{m  ), 

generic-l_surface_area(m  ),generic-2-surface_area-(m  ), material-code,  density _(kg/m  ), 

specific_heat_(J/kgC),  absorptivity,  emissivity,  thermal-conductivity _(W/mK),  all-adjacent-regions 


3.3.  Firpass  Examples 


3.3.1.  PRISM  File  Example 
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In  the  following  examples  the  mged  file  test.veh.g  will  be  used  as  input  to  firpass.  This  file  pertains  to  the 
model  of  a  simple  vehicle  developed  for  use  in  the  following  examples.  The  following  is  an  example  where 
the  facet  records  of  a  PRISM  file  are  created  using  this  mged  file.  Program  prompts  and  messages  are 
printed  in  italics.  User  response  is  printed  in  an  upright  font. 


$  firpass  test.veh.g  vehicle  air 

Write  output  to  standard  out  (0)  or  a  file  (1)  or  not  at  all  (2)9  1 
Enter  name  of  file  output  is  to  be  written  to  (15  char  max),  veh.f.out 
Enter  name  of  error  file  to  be  created  (15  char  max),  veh.f.err 
Enter  name  of  second  pass  file  to  be  created  (15  char  max).  veh.f.2nd 
Enter  name  of  region  ^  ^  name  file  to  be  created  (15  char  max),  veh.f.rnn 
Enter  name  of  mater-al  id  file  to  be  read  (15  char  max),  prp 
Enter  type  of  file  to  be  written. 

0  -  PRISM  file 

1  -  Generic  file 

2  -  Geometric  properties  file 

0 

Enter  name  of  facet  file  to  be  created.  (15  char  max)  veh.f.fac 
Which  release  of  PRISM  is  being  used,  2.0  (2)  or  S.O  (S)9  3 
File  Used:  test.veh.g 
Material  ID  File:  prp 

Database  Title:  test  vehicle  for  use  with  irprep  programs 
Mallocing  arrays. 

Zeroing  variables. 

Center  of  bounding  rpp  (356.500000,  0.000000,  1250.000000  ) 

Length  of  diagonal  of  bounding  rpp:  3829.051456 
Mtnimums  &  maximums 
x:  -3472.551456  -  4185.551456 
y:  -3829.051456  -  3829.051456 
z:  -2579.051456  -  5079.051456 
Model  minimum  &  maximum. 

X:  -3001.000000  to  3714.000000 
Y:  -1350.000000  to  1350.000000 
Z:  0.000000  to  2500.000000 

Enter  grid  spacing  (mm)  for  fired  rays. 

10 

SHOOTING  DOWN  X-AXIS 

SUMMARY  OF  FILES  USED  &  CREATED 
.g  file  used:  test.veh.g 
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regions  used: 

vehicle 

air 

material  id  file  used:  prp 
output  file  created-  veh.f.out 
second  pass  file  created:  veh.f.Snd 
error  file  created:  veh.f.err 
region  #  &  name  file  created:  veh.f.rnn 
facet  file  created:  veh.f.fac 
(format  is  PRISM  3.0) 


Freeing  memory. 

$ 


The  command  line  for  using  firpass  is 
firpass  model. g  objects... 

where  model. g  is  the  mged  geometry  to  be  used  (test.veh.g)  and  objects  are  all  the  groups  to  be  interrogated 
(vehicle  and  air).  Note  that  for  this  particular  run  the  objects  vehicle  and  air  were  used  in  that  order.  This 
will  allow  the  air  to  be  stripped  out  easily  where  it  is  not  needed.  The  objects  are  ray-traced  in  the  order 
they  are  given;  therefore,  air  will  appear  at  the  end  of  the  file  and  less  renumbering  of  regions  will  need  to  be 
done.  For  this  particular  run  the  user  chose  to  create  an  output  file,  naming  it  veh.f.out.  The  .f.out  suffix 
means  that  it  is  the  output  file  created  by  firpass.  This  convention  for  naming  files  will  be  used  throughout 
the  manual.  The  user  also  needed  to  enter  the  name  of  an  error  file  (.f.err  implies  the  error  file  created  by 
firpass).  The  user  then  entered  veh.f.2nd  for  the  second  pass  file.  This  file  will  be  used  with  secpass,  while 
veh.f.rnn  will  be  used  with  shapefact  and  display.  The  type  of  file  to  be  created  is  the  "facet"  records  for 
PRISM;  therefore,  the  user  chooses  0  (PRISM  file).  A  facet  file  will  now  be  created.  The  version  of  PRISM 
to  be  used  must  be  chosen.  In  this  example  version  3.0  was  selected.  The  formatting  between  the  2.0  and 
3.0  version  of  PRISM  is  slightly  different.  The  name  of  a  material  file  needed  to  be  entered,  prp.  In  this 
particular  material  properties  file  absorptivity  was  set  to  0.60  and  emissivity  was  set  to  0.94.  The  name 
(including  path  information)  must  be  less  than  or  equal  to  fifteen  characters.  Information  about  the  file  and 
the  size  of  the  geometry  is  printed  out  in  millimeters  and  the  user  must  choose  a  grid  spacing.  The 
geometry  example  is  approximately  six  meters  long,  two  meters  wide,  and  two  and  a  half  meters  high.  A 
grid  spacing  of  10  millimeters  was  chosen  as  being  an  appropriate  spacing.  Since  the  run  took  less  than  10 
minutes  this  spacing  could  be  made  smaller  and  still  run  in  a  reasonable  period  of  time.  The  program  prints 
SHOOTING  DOWN  X-AXIS  as  rays  are  fired  even  though  the  rays  are  actually  being  fired  in  an  arbitrary 
direction.  Firpass  ends  with  a  summary  of  the  files  used  and  created.  A  copy  of  veh.f.fac  created  by 
running  this  sample  case  may  be  found  in  Appendix  B. 

Let’s  look  at  one  region  in  the  mged  file  and  examine  the  values  computed.  Region  r.hull03  is  the  back  wall 
of  the  hull,  it  is  bounded  on  the  rear  by  exterior  air  (nothing)  and  on  the  front  by  engine  air  (air  type  5). 


The  region  number  of  r.huilOS  is  three,  this  can  be  found  in  file  veh.f.fac.  This  region  is  rectangular  and 
ranges  from  -3000  to  -2980  in  the  x-direction,  -980  to  980  in  the  y-direction,  and  550  to  1650  in  the  z- 
direction  (units  are  millimeters).  The  material  id  is  1.  With  this  information  the  exterior  surface  area  may 
be  computed  by  hand  (2.156m  )  and  the  mass  may  be  computed  by  hand  (335kg).  The  surface  unit  normal 
should  be  (-1,  0,  0).  These  values  match  very  v/ell  with  the  values  found  by  firpass  as  seen  in  veh.f.fac. 

3.3.2.  Generic  File  Example 

The  following  is  a  run  of  firpass  where  a  generic  file  is  created  using  the  same  vehicle  as  before. 

$  firpass  test.veh.g  vehicle  air 

Write  output  to  standard  out  (0)  or  a  file  (l)  or  not  at  all  (2)9  2 

Enter  name  of  error  file  to  be  created  (15  char  max),  veh.f.err 

Enter  name  of  region  #  &  name  file  to  be  created  (15  char  max),  veh.f.rnn 

Enter  name  of  second  pass  file  to  be  created  (15  char  max).  veh.f.2nd 

Enter  name  of  material  id  file  to  be  read  (15  char  max),  prp 

.Cnter  type  of  file  to  be  written. 

0  ■  PRISM  file 

1  -  Generic  file 

2  -  Geometric  properties  file 

1 

Enter  name  of  generic  file  to  be  created.  (15  char  max)  veh.f.gen 
File  Used:  test.veh.g 
Material  ID  File:  prp 

Database  Title:  test  vehicle  for  use  with  irprep  programs 
Mallocing  arrays. 

Zeroing  variables. 

Center  of  bounding  rpp  (  S56. 500000,  0.000000,  1250.000000  ) 

Length  of  diagonal  of  bounding  rpp:  3829.051456 
Minimums  &  maximums 
x:  -3472.551456  -  4185.551456 
y:  -3829.051456-  3829.051456 
z:  -2579.051456-  5079.051456 
Model  minimum  &  maximum. 

X:  -3001 . 000000  to  3714. 000000 
Y:  -1350.000000  to  1350.000000 
Z:  0.000000  to  2500.000000 

Enter  grid  spacing  (mm)  for  fired  rays. 
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SHOOTING  DOWN  X-AXIS 

10 


SUMMARY  OF  FILES  USED  &  CREATED 
.g  file  used:  test.veh.g 
regions  used: 

vehicle 

air 

material  id  file  used:  prp 

second  pass  file  created:  veh.f.Snd 

error  file  created:  veh.f.err 

region  #  &  name  file  created:  veh.f.rnn 

generic  file  created:  veh.J.gen 


Freeing  memory. 

$ 


The  only  difference  between  this  and  the  previous  example  is  that  the  user  chose  to  create  a  generic  file.  The 
generic  file,  veh.f.gen,  may  be  found  in  Appendix  B. 


3.3.3.  Geometric  Properties  File  Example 

Below  is  a  sample  session  where  the  geometric  properties  file  is  to  be  written  by  firpass,  using  the  same 
geometry  as  in  the  last  two  sections. 


$  firpass  test.veh.g  vehicle  air 

Write  output  to  standard  out  (0)  or  a  file  (l)  or  not  at  all  (2)9  2 
Enter  name  of  error  file  to  be  created  (15  char  max),  veh.f.err 
Enter  name  of  second  pass  file  to  be  created  (15  char  max).  veh.f.2nd 
Enter  name  of  region  #  &  name  file  to  be  created  (15  char  max),  veh.f.rnn 
Enter  name  of  material  id  file  to  be  read  (15  char  max),  prp 
Enter  type  of  file  to  be  written. 

0-  PRISM  file 

1  -  Generic  file 

2  -  Geometric  properties  file 

2 

Do  you  want  a  readable  (0)  or  non-readable  (1)  geometric  file?  0 
Enter  name  of  geometric  properties  file  to  be  created  (15  char  max),  veh.f.geo 
File  Used:  test.veh.g 
Material  ID  File:  prp 

Database  Title:  test  vehicle  for  use  with  irprep  programs 
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Mallocing  arrays. 

Zeroing  variables. 

Center  of  bounding  rpp  (  S56. 500000,  0.000000,  1S50. 000000  ) 
Length  of  diagonal  of  bounding  rpp:  3829.051456 
Minimums  &  maximums 
x:  -3472.551456  ■  4185.551456 
y:  -3829.051456  -  3829.051456 
z:  -2579.051456  -  5079.051456 
Model  minimum  &  maximum. 

X:  -3001.000000  to  3714.000000 
Y:  -1350.000000  to  1350.000000 
Z:  0.000000  to  2500.000000 

Enter  grid  spacing  (mm)  for  fired  rays 
10 

SHOOTING  DO WN X-AXIS 

SUMMARY  OF  FILES  USED  &  CREATED 
.g  file  used:  lest.veh.g 
regions  used: 

vehicle 

air 

material  id  file  used:  prp 

second  pass  file  created:  veh.f.2nd 

error  file  created:  veh.f.err 

region  #  &  name  file  created:  veh.f.rnn 

geometric  file  created:  veh.f.geo 


Freeing  memory. 

$ 

The  geometric  properties  file,  veh.f.geo,  may  also  be  found  in  Appendix  B. 


4  SECPASS 

4.1.  Secpass  Input  Files 

Secpass  interrogates  an  mged  file  to  obtain  the  conductivity  table  (conduction  factors  for  adjacent  regions) 
for  a  PRISM-ready  file  or  a  generic  file.  When  the  "facet"  records  from  firpass  are  combined  with  the 
conductivity  table  of  secpass,  a  PRISM-ready  file  can  be  made.  A  compiled  copy  of  secpass,  the  same 
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mged  file  used  with  firpass,  the  same  material  properties  file  used  with  firpass,  and  the  second  pass  file 
created  by  firpass  for  use  with  secpass  are  needed  to  run  secpass.  A  second  pass  file  from  firpass  is 
generated  regardless  of  which  output  file  one  selects.  The  second  pass  file  contains  the  total  number  of 
regions  in  the  model.  An  error  will  occur  if  this  number  of  regions  does  not  equal  the  number  of  regions  in 
the  mged  geometry  file  named  in  the  command  line.  For  each  region  listed  in  the  second  pass  file  there  will 
be  a  region  number,  centroid,  material  id,  and  shared  surface  area  for  each  adjacent  region. 


4.2  Secpass  Output  Files 

There  will  be  three  or  four  output  files  generated  by  secpass;  an  optional  diagnostic  file,  an  error  file,  an 
intermediate  file,  and  a  primary  file.  The  diagnostic  file  is  optional,  used  mainly  for  debugging  purposes.  It 
may  be  written  to  a  file,  written  to  the  screen,  or  not  written  at  all.  The  error  file  will  contain  information 
on  any  adjacent  regions  for  which  the  number  of  calculations  (number  of  rays  that  hit  these  two  regions)  to 
determine  the  conductance  is  less  then  the  minimum  number  set  in  secpass. c.  The  intermediate  file 
contains  the  four  conductance  values  between  each  adjacent  region  calculated  using  four  types  of  conduction 
path  length  approximations.  The  four  conduction  path  length  approximations  are:  average,  root  mean 
squared  (rms),  minimum,  and  maximum.  The  user  may  compare  the  conductance  values  between  adjacent 
regions  for  different  conduction  path  lengths  and  determine  which  is  the  most  appropriate. 

The  user  creates  one  primary  file:  a  conductance  file  or  generic  file.  The  conductance  file  contains  the 
information  needed  by  PRISM.  The  format  of  this  file  may  be  found  in  the  PRISM  User’s  Manual  [3,4]. 
The  information  contained  in  this  file  for  use  with  PRISM  is  the  region  number,  the  adjacent  region 
number,  the  conductance  between  the  two  regions,  and  the  shared  surface  area.  The  shared  surface  area  is 
not  actually  needed  by  PRISM  but  is  printed  so  that  the  user  may  see  what  the  area  is.  This  file  is  written 
out  in  a  triangular  format,  herein  the  conductance  between  regions  m  and  n  is  not  printed  again  as  a 
conductance  between  regions  n  and  m. 

The  user  may  choose  to  create  a  generic  file  instead  of  the  conductance  file.  The  generic  file  option  is  used  in 
conjunction  with  other  thermal  models.  The  user  must  extract  and  format  the  information  needed  by  his 
thermal  model.  For  each  region  in  the  mged  geometry  model  the  following  information  will  be  printed  in 
the  generic  file. 

4  region-number  number_of_adjacent_regions 

2 

adjacent-region-number  shared_area-(m  )  conduction_distance_(m) 

The  number  four  indicates  that  the  following  lines  will  contain  information  on  conduction  between  two 
adjacent  regions.  The  number  of  adjacent  regions  defines  how  many  lines  of  information  follow  the  current 
line. 


4.3.  Secpass  Examples 


4.3.1  PRISM  File  Example 
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The  following  example  uses  the  same  mged  model  as  before. 


$  secpass  test.veh.g  vehicle  air 

Write  output  to  standard  out  (0)  or  a  filefl)  not  at  all  (S)9  1 
Enter  name  of  output  file  (15  char  max),  veh.s.out 
Enter  name  of  file  that  has  second  pass  information 
in  it  (15  char  max).  veh.f.2nd 
Enter  type  of  output  file  to  be  generated. 

0  ■  PRISM  File 
1  -  Generic  File 

0 

Enter  name  of  file  to  be  created  for  PRISM  conductivity 

information  (15  char  max),  veh.s.con 

Which  release  of  PRISM  is  being  used,  2.0  (2)  or  S.O  (S)9  3 

Which  length  calculation  should  be  used  when 

computing  conduction 

between  regions? 

1  -  average  length 

2  -  rms  length 

S  -  minimum  length 
-  maximum  length 
3 

Enter  name  of  file  to  be  created  for  conductivity 

table  (15  char  max),  veh.s.tbl 

Enter  name  of  material  file  (15  char  max),  prp 

Enter  name  of  error  file  to  be  created  (15  char  max),  veh.s.err 

Database  title:  test  vehicle  for  use  with  irprep  programs 

Number  of  regions  in  mged  file:  21 

Mallocing  arrays. 

cond  malloced 

loop  malloced 

All  variables  zeroed. 

second  pass  file  opened 

The  number  of  regions  read  was  21 

Center  of  bounding  rpp  (  S56. 500000,  0.000000,  1250.000000  ) 
Length  of  diagonal  of  bounding  rpp:  3829.051456 
Minimums  &  maximums  of  grid: 

■3472.551456  ■  4185.551456 
-3829.051456  -  3829.051456 
■2579.051456  ■  5079.051456 

Model  minimum  &  maximum. 

X:  -3001.000000  to  3714.000000 
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Y:  -1350.000000  to  1350.000000 
Z:  0.000000  to  2500.000000 


Enter  spacing  (mm)  between  fired  rays.  10 
grid  spacing:  10.000000 
Shooting  down  x-axis. 


SUMMARY  OF  FILES  USED  &  CREATED 
.g  file  used:  test.veh.g 
regions  used: 

vehicle 

air 

file  containing  second  pass  information:  vch.f.2nd 
material  file  used:  prp 
output  file  created:  veh.s.out 
conductivity  file  created:  veh.s.con 
(format  is  PRISM  3.0) 
conductivity  table  file  created:  veh.s.tbl 
error  file  created:  veh.s.err 


Freeing  memory. 

% 


The  command  line  for  using  seepass  is 
seepass  model. g  objects... 

where  model. g  is  the  mged  geometry  file  to  be  used,  and  objects  are  all  groups  to  be  interrogated.  The  user 
must  list  the  objects  in  the  same  order  as  the  objects  were  listed  when  using  firpass.  If  the  objects  are  not 
listed  in  the  same  order,  the  region  numbers  will  not  correspond  to  the  region  numbers  in  the  files  created 
with  firpaas.  For  this  particular  run  the  user  chose  to  create  an  output  file,  naming  it  veh.s.out  (.s.out 
implies  the  output  file  created  by  seepass).  The  user  now  decides  to  create  a  conductance  PRISM  file  and 
must  enter  a  name  for  the  file,  veh.s.con  (  s  con  implies  the  conductance  file  for  PRISM  created  by  seepass). 
The  user  is  now  asked  which  version  of  PRISM  is  to  be  used,  2.0  or  3.0.  There  is  a  slight  difference  in  the 
format  of  these  two  files  [3,4].  In  this  particular  example  the  user  is  using  PRISM  3.0.  The  user  selects  the 
conduction  path  length  form  to  be  used  when  computing  the  conductance  between  regions.  The  average 
length  is  the  average  length  of  ray  segments  traced  between  the  centroid  and  shared  surface.  The  rms 
length  is  the  root  mean  squared  length  of  these  ray  segments  between  the  centroid  and  shared  surface,  and 
the  minimum  and  maximum  lengths  are  the  minimum  and  maximum  lengths  of  the  ray  segments.  In  this 
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particular  example  the  minimum  length  was  chosen.  The  user  now  enters  a  file  name  for  the  intermediate 
file,  veh.s.tbl  (.s.tbl  implies  intermediate  conductance  table  created  by  secpass).  This  table  will  contain  all 
four  conduction  path  lengths  and  conductances  between  adjacent  regions  for  all  regions  in  the  geometry  file 
to  assist  the  user  in  choosing  a  conduction  path  length  form.  Next  the  user  enters  the  material  properties 
file,  the  same  file  used  with  firpass.  Again,  an  error  file  is  created  containing  information  on  regions  where 
the  conduction  path  length  was  calculated  using  fewer  than  a  minimum  number  of  values.  Currently  this 
minimum  number  of  values  (variable  MINCAL)  in  secpass.c  is  set  to  ten  which  represents  the  number  of 
rays  intersecting  the  shared  surface.  Next,  information  about  the  database  is  given  such  as  the  title  and 
number  of  regions.  The  next  five  statements  simply  tell  where  the  program  is  in  its  execution  and  are 
probably  of  little  interest  to  the  user.  The  next  line  indicates  the  number  of  regions  in  the  file  created  by 
firpass.  If  this  number  does  not  agree  with  the  number  of  regions  in  the  BRL-CAD  geometry  file,  an  error 
will  occur.  The  next  lines  are  information  on  the  dimensions  of  the  model  in  millimeters.  The  user  needs  to 
enter  the  grid  spacing  for  the  rays  that  are  to  be  fired.  This  spacing  should  be  the  same  as  that  used  in 
firpass.  When  the  run  is  finished,  a  summary  of  the  files  used  and  created  will  be  printed.  The  file 
veh.s.con  may  be  found  in  Appendix  C. 

The  conductance  between  two  regions  may  be  checked  by  hand.  Regions  r.hull03  (region  number  3)  and 
r.hull04  (region  number  4)  are  next  to  each  other.  They  are  both  rectangular  regions.  r.hull04  is  the 
bottom  plate  of  the  hull  and  r.hull03  is  the  back  plate  of  the  hull.  The  x-component  of  r.hull04  extends 
from  -3000  to  3000,  the  y-component  extends  from  -980  to  980,  and  the  z-component  extends  from  550  to 
570  (millimeters).  The  x-component  of  r.hull03  extends  from  -3000  to  -2980,  the  y-component  extends  from 
-980  to  980,  and  the  z-component  extends  from  570  to  1630.  The  centroid  of  r.hullOS  is  (-2990,  0,  1100)  and 
the  centroid  of  r.hull04  is  (0,  0,  560).  Both  these  centroids  agree  with  what  was  found  in  firpass.  The 
minimum  distance  from  the  centroid  of  r.hull04  to  the  shared  surface  area  is  the  distance  from  the  centroid 
to  where  the  right  shared  surface  area  of  r.hull03  meets  r.hull04  (2.98m).  The  minimum  distance  from  the 
centroid  of  r.hull03  to  the  edge  of  r.hull04  is  the  distance  straight  down  fro  the  centroid  to  the  top  of 
r.hull04  (0.53m).  Both  of  these  agree  with  what  was  found  using  secpass.  The  shared  surface  area  is  0.04 

O 

m"  Both  regions  are  mild  steel  and  have  a  thermal  conductivity  of  52.000  W/mK.  Now  the  equation  from 
the  FRI.SM  manual  may  be  used  (3,4). 

_  f  J:l  iL  ” 

Ik, A  K,A 

RK;  j  :  conductance  between  regions  i  and  j 
Lj  :  conduction  path  length  of  region  i 
Lj  :  conduction  path  length  of  region  j 
Kj  :  thermal  conductivity  of  region  i 
Kj  :  thermal  conductivity  of  region  j 
A:  shared  surface  area  between  region  i  and  j 

When  these  values  are  substituted  in  the  equation  above,  the  conductance  between  regions  r.hull03  and 
r.hull04  is  found  to  be  0.593  which  agrees  with  the  value  calculated  by  secpass. 
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4.3.2.  Generic  File  Example 


The  following  is  the  same  run  as  in  the  last  section  except  that  a  generic  file  was  created  instead  of  a  PRISM 
file. 


$  secpass  test.veh.g  vehicle  air 

Write,  output  to  standard  out  (0)  or  a  file(l)  not  at  all  (2)9  2 
Enter  name  of  file  that  has  second  pass  information 
in  it  (15  char  max).  veh.f.2nd 
Enter  type  of  output  file  to  be  generated. 

0  -  PRISM  File 
1  -  Generic  File 

1 

Enter  name  of  generic  file  to  be  created  (15  char  max),  veh.s.gen 
Which  length  calculation  should  be  used  when 
computing  conduction 
between  regions? 

1  -  average  length 

2  -  rms  length 

S  -  minimum  length 
4  •  maximum  length 
3 

Enter  name  of  file  to  be  created  for  conductivity 

table  (15  char  max),  veh.s.tbl 

Enter  name  of  material  file  (15  char  max),  prp 

Enter  name  of  error  file  to  be  created  (15  char  max),  veh.s.err 

Database  title:  tied  MGED  Database 

Number  of  regions  in  mged  file:  21 

Mallocing  arrays. 

cond  malloced 

loop  malloced 

All  variables  zeroed. 

second  pass  file  opened 

The  number  of  regions  read  was  21 

Center  of  bounding  rpp  (  S56. 500000,  0.000000,  1250.000000  ) 
Length  of  diagonal  of  bounding  rpp:  S829. 051456 
Minimums  &  maximums  of  grid: 

■3472.551456  -  4185.551456 
-3829.051456  -  3829.051456 
■2579.051456  ■  5079.051456 

Model  minimum  &  maximum. 

X:  -3001.000000  10  3714.000000 
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Y:  -1S50.000000  to  1S50. 000000 
Z:  0.000000  to  2500.000000 


Enter  spacing  (mm)  between  fired  rays.  10 

grid  spacing:  10.000000 

Shooting  down  i-axis. 

Opened  generic  file. 


SUMMARY  OF  FILES  USED  &  CREATED 

.g  file  used:  test.veh.y 
regions  used: 

vehicle 

air 

file  containing  second  pass  information:  veh.f.2nd 
material  file  used:  prp 
generic  file  created:  veh.s.gen 
conductivity  table  file  created:  veh.s.tbl 
error  file  created:  veh.s.err 


Freeing  memory. 
$ 


The  user  chose  to  create  a  generic  file  this  time.  All  the  other  answers  to  the  questions  are  as  before.  The 
file  veh.s.gen  may  be  found  in  Appendix  C. 


5  SHAPEFACT 

0.1.  Shapefact  Input  Files 

Shapefact  interrogates  an  mged  file  to  find  the  shape  (intervisibility)  factors  between  regions,  primarily  for 
.Jim  Rapp’s  engine  compartment  radiation  models  [5).  In  order  to  run  shapefact  the  user  will  need  a 
compiled  copy  of  shapefact,  an  mged  geometry  file  of  engine  components,  and  a  region  number  and  name 
file  created  by  firpass.  This  geometry  file  must  have  engine  air  correctly  placed  and  labeled.  Shapefact 
only  uses  rays  for  finding  the  shape  factor  if  they  leave  a  region  and  enter  engine  air.  In  general  shapefact 
finds  the  shape  factors  of  an  engine  which  is  a  subset  of  a  larger  mged  geometry  file.  Usually  the  mged 
geometry  file  used  with  firpass  and  secpass  is  an  entire  vehicle  including  the  engine.  Therefore,  the  region 
number  and  name  file  created  by  firpass  is  needed  so  that  the  regions  of  the  engine  being  interrogated  by 
shapefact  have  the  same  region  iiumbers  as  the  engine  regions  in  the  original  model. 
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5.2.  Shapefact  Output  Files 


There  will  only  be  one  of  two  types  of  primary  files  and  an  error  file  generated:  a  regular  or  generic  file. 
The  regular  file  is  for  use  with  Jim  Rapp’s  engine  compartment  radiation  models  [5].  In  this  file  the  first 
line  contains  the  number  of  forward  rays  fired.  The  succeeding  lines  will  contain  information  about  each 
region.  The  format  is  below. 

region-number  region-name  engine-area-(m  ) 
region-number  shape-factor 
sum— of-hits 
sum— of-shape_f  actors 

The  sum  of  hits  and  the  sum  of  shape  factors  is  provided  for  error  checking.  For  an  enclosed  engine  the  sum 
of  the  shape  factors  shoul  1  sum  to  1.0.  The  generic  file  has  been  ’vritten  to  be  used  with  other  engine 
models.  It’s  format  is  shown  below. 

2 

5  region-number  radiation-surface_area_(m  }  number-of-viewable-regions 
viewable-region-number  shape-factcr 

The  5  indicates  the  data  will  be  radiat’'>n  i  ^at  transfer  dat-’  I’hp  iVamber  of  viewable  regions  indicates  how 
many  regions  the  current  reg'on  has  shape  •“'■.tors  with,  this  also  indicates  the  number  of  data  lines  to 
follow. 


5.3.  Shapefact  Examples 
5  3.1.  a, pie  One:  Sphere  and  Disk 

The  geometry  file  in  this  example  consists  of  a  sphere  and  disk  enclosed  in  a  hollow  box.  Any  empty  space 
will  be  filled  with  engine  air. 

Below  is  the  script  for  running  this  particular  example. 


$  shapefact  sph.dsk.g  all. air 

Enter  type  of  file  to  be  written  (0=  ^regular  or  1=  ^  generic).  0 
Enter  name  of  ovlput  file  (15  char  max). 
sd.sh.reg 

Enter  the  name  of  the  error  file  (15  char  max). 
sd.sh.err 

Enter  region  ff  &  name  file  to  be  read  (15  char  max). 
sd.f.rnn 

Do  you  want  to  dump  tnlerm  shape  factors  to  screen  (0-no,  I-yes)?  0 
Enter  number  of  rays  to  be  fired.  5000000 
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Do  you  wish  to  enter  your  own  seed  (0)  or  use  the  default  of  1  (l)9  1 
seed  initialized 

Database  Title:  sphere  &  disk  for  use  with  shape 
Number  of  regions:  4 

1000000.000000  rays  have  been  fired  in  forward  direction. 
2000000.000000  rays  have  been  fired  in  forward  direction. 

SOOOOOO. 000000  rays  have  been  fired  in  forward  direction. 
4000000.000000  rays  have  been  fired  in  forward  direction. 
5000000.000000  rays  have  been  fired  in  forward  direction. 

Finding  correct  region  numbers. 

Finished  finding  correct  region  numbers. 

$ 


Shapefact  is  started  by  typing  the  following  line, 
shapefact  model. g  objects... 

The  first  information  requested  is  the  type  of  file  to  be  created  and  the  name  of  the  output  file  where  all  the 
information  on  shape  factors  will  be  recorded  as  well  as  the  name  of  an  error  file  to  be  created.  The  user 
also  needs  to  enter  a  region  number  and  name  file  that  was  created  by  firpass.  This  is  done  so  the  region 
numbering  of  the  engine  corresponds  to  the  numbering  of  the  entire  vehicle  file.  If  the  user  so  chooses, 
intermediate  shape  factors  may  be  printed  out.  If  the  user  chooses  to  print  out  intermediate  shape  factors, 
they  will  be  printed  every  million  rays.  Columns  one  and  two  contain  region  numbers  and  column  three 
contains  the  shape  factor  between  the  two  regions.  Next  the  user  needs  to  enter  the  number  of  rays  to  be 
fired.  This  will  be  a  VERY  large  number  since  the  majority  of  the  rays  will  miss  the  model  altogether. 
Then  the  user  may  choose  to  enter  his  own  seed  for  the  random  number  generator.  The  regular  file  sd.sh.reg 
may  be  found  in  Appendix  D.  The  questions  and  answers  for  creating  a  generic  file  are  exactly  the  same 
except  the  user  must  provide  an  answer  of  1  for  the  first  question.  A  generic  file  sd.sh.gen  may  also  be  found 
in  Appendix  D. 

The  sphere  in  this  example  has  a  radirs  of  500mm  and  is  centered  at  (0,  0,  500).  The  disk  has  a  radius  of 
1000mm  and  the  top  face  is  centered  at  (0,  0,  -500).  The  shape  factor  for  this  geometry  may  be  computed 
using  the  following  equation  [8]. 


2 

“I  \ 

1 

'rd  ' 
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1  - 

1  + 
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.  h  > 

:  shape  factor  from  sphere  to  disk 
rj  :  radius  of  disk 

h:  distance  from  center  of  sphere  to  top  face  of  disk 


Hand  computation  of  the  shape  factor  from  the  sphere  to  the  disk  yields  0.146,  which  agrees  with  the  output 
from  shapefact.  From  this  output,  one  infers  that  the  shape  factor  from  region  number  two  (r.sph,  sphere) 
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to  region  number  one  (r.disk,  disk)  is  0.146196. 


5.3.2.  Example  Two:  Concentric  Spheres 

This  example  finds  the  shape  factors  for  two  concentric  spheres  enclosed  in  a  sphere.  The  two  concentric 
spheres  are  enclosed  in  a  sphere  so  that  rays  leaving  the  spheres  are  not  taken  into  account. 

$  shapefact  con.sph.g  all. air 

Enter  type  of  file  to  be  written  ^regular  or  l=>generic).  0 
Enter  name  oj  output  file  (15  char  max). 
cs.sh.reg 

Enter  the  name  of  the  error  file  (15  char  max). 
cs.sh.err 

Enter  the  region  #  &  name  file  to  be  read  (15  char  max). 
cs.f.rnn 

Do  you  want  to  dump  interm  shape  factors  to  screen  (0-no,  1-yes)?  0 
Enter  number  of  rays  to  be  fired.  5000000 

Do  you  wish  to  enter  your  own  seed  (0)  or  use  the  default  of  1  (l)?  1 
seed  initialized 

Database  Title:  concentric  spheres  for  use  with 
Number  of  regions:  4 

1000000.000000  rays  have  been  fired  in  forward  direction. 

2000000.000000  rays  have  been  fired  in  forward  direction. 

3000000.000000  rays  have  been  fired  in  forward  direction. 

4000000 .000000  rays  have  been  fired  in  forward  direction. 

5000000 .000000  rays  have  been  fired  in  forward  direction. 

Finding  correct  region  numbers. 

Finished  finding  correct  region  numbers. 

$ 

The  regular  file,  cs.sh.reg,  may  be  found  in  Appendix  D  along  with  cs.sh.gen,  a  generic  file.  Again,  this 
calculation  may  be  checked  by  hand  using  the  following  equation  jS]. 

rj  :  radius  of  disk.  257text:  S12  <- 12;  b=0,h=42,lf=l,rf=l 


F,2  :  shape  factor  from  sphere  one  to  sphere  two 

Fji  :  shape  factor  from  sphere  two  to  sphere  one 

V22  '■  shape  factor  from  sphere  two  to  sphere  two 
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r,  :  radius  of  sphere  one 
rj  :  radius  of  sphere  two 


Sphere  one  is  the  inside  sphere  and  sphere  two  is  the  outside  sphere.  The  radius  of  sphere  one  is  1000  mm 
and  of  sphere  two  is  2000  mm.  The  shape  factor  from  sphere  one  to  sphere  two  is  1.0,  the  shape  factor  from 
sphere  two  to  sphere  one  is  0.25,  and  the  shape  factor  from  sphere  two  to  itself  is  0.75.  These  values  agree 
fairly  well  with  shapefact.  The  output  from  shapefact  states  that  the  shape  factor  from  sphere  one 
(region  1,  r.sphl)  to  sphere  two  (region  2,  r.sph2)  is  one.  The  shape  factor  from  sphere  two  to  sphere  one  is 
0.249554.  And  the  shape  factor  from  sphere  two  to  itself  is  0.750446.  If  closer  agreement  is  needed,  more 
rays  may  be  fired.  However,  since  the  method  of  calculation  involves  Monte  Carlo  sampling,  to  increase  the 
accuracy  of  the  current  answer  by  one  more  significant  digit  required  100  times  more  ray-tracing! 


6.  CREATING  A  PRISM  VEHICLE  FILE  FROM  IRPREP  OUTPUT  FILES 
6.1.  Modifications  to  Irprep  Output  Files 

A  PRISM  vehicle  file  is  created  by  combining  the  facet  file  created  by  firpass  and  the  conductivity  file 
created  by  secpass  into  a  vehicle  file.  PRISM  does  not  accept  air  regions;  therefore,  all  air  regions  must  be 
removed  from  the  facet  records  and  any  reference  to  them  must  be  removed  from  the  conductivity  table. 
The  file  must  be  renumbered  so  the  region  numbers  are  consecutive.  This  is  why  it  is  recommended  that  the 
air  be  a  separate  group  placed  last  when  running  firpass  and  secpass.  At  the  very  least  the  last  line  of  the 
facet  records,  indicating  the  end  must  be  renumbered.  The  first  line  of  the  conductivity  table  indicating  this 
file  is  for  use  with  PRISM  must  be  removed.  Firpass  labels  all  region  types  one.  Any  special  regions  such 
as  tracks  and  wheels  need  to  be  relabled.  (See  references  3  and  4  for  a  complete  listing  of  special  regions.) 
Any  track  regions  must  be  labeled  right  or  left,  also  [3,4].  And  wheels  need  to  have  a  radius  and  if  the 
vehicle  is  moving  a  friction  associated  with  them.  Any  regions  that  were  not  divided  correctly  should  have 
their  normals  altered.  For  example,  if  the  gun  tube  was  one  region  the  normal  should  probably  be  pointing 
straight  up,  (0,0,1).  Any  wheels  that  are  one  region  should  probably  have  a  normal  of  (0,1,0)  or  (0,-l,0) 
depending  whether  they  are  on  the  left  or  right. 


6.2.  An  Example 

First  veh.f.fac  and  veh.s.con  are  combined  to  create  a  PRISM  vehicle  file,  veh.veh.  Remove  the  first  line  of 
what  was  veh.s.con,  "Conductivity  file  for  use  with  PRISM.”  This  line  was  just  printed  to  let  the  user  know 
how  the  file  was  to  be  used. 

Next  remove  all  lines  in  the  facet  records  that  deal  with  air.  For  this  example  that  is  regions  nineteen 
through  twenty-one.  The  last  line  of  the  facet  records  needs  to  be  renumbered  nineteen.  In  the  conductivity 
table  part  of  the  file  all  references  to  regions  nineteen  through  twenty-one  (air  regions)  need  to  be  removed. 

Now  any  special  region  types  need  to  be  identified.  The  only  special  regions  in  this  file  are  the  wheels, 
regions  fourteen  through  seventeen.  These  regions  are  of  type  eleven;  therefore,  type  one  must  be  changed 
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to  eleven.  Regions  twelve  and  thirteen  are  part  of  the  group  wheels  but  they  are  axles  not  wheels;  therefore, 
the  type  does  not  need  to  be  changed.  When  a  region  is  identified  as  a  wheel  it  must  have  an  associated 
radius.  The  radius  variable  for  this  example  needs  to  be  set  to  0.5  since  the  wheel  radius  is  a  half  meter. 
The  friction  variable  was  not  changed  since  in  this  scenario  the  vehicle  will  not  be  moving. 

Any  normals  that  need  to  be  adjusted  may  be  changed  now.  Since  the  wheels  are  all  one  region  the  left 
wheels  will  get  a  normal  of  (0,1,0)  and  the  right  wheels  will  get  a  normal  of  (0,-l,0).  The  guntube  was  also 
modeled  as  one  region  so  the  normal  was  chosen  to  be  (0,0,1). 

The  vehicle  file  is  now  ready  to  be  run  through  PRISM  with  a  set  of  files  that  were  supplied  with  PRISM: 
sc071984,  statdrv,  w071984,  and  tr071984.  These  are  a  scenario  file,  a  stationary  driver  file,  a  weather  file, 
and  a  terrain  file. 

A  temperature  output  file,  veh.prm,  was  created.  This  file  contains  the  temperature  of  each  region  at 
different  time  steps.  See  the  PRISM  manual  [3,4]  for  the  format  of  this  file.  The  vehicle  file,  veh.veh;  the 
temperature  output  file,  veh.prm;  and  the  PRISM  input  files  may  be  found  in  Appendix  E. 


7  DISPLAY 

7.1.  Display  Input  Files 

Display  ray-traces  a  BRL-CAD  geometry  file.  As  a  ray  hits  a  region,  the  temperature  of  that  region  is 
looked  up  in  the  PRISM  file  and  a  temperature  is  assigned  to  the  pixel  corresponding  to  that  ray.  In  order 
to  run  display  a  compiled  copy  of  display,  the  BRL-CAD  geometry  file,  the  temperature  output  file  from 
PRISM  (or  a  generic  temperature  output  file),  and  the  region  number  and  name  file  created  by  firpass  must 
be  available.  In  addition  to  this  the  user  must  know  the  number  of  regions  (if  using  a  PRISM  output  file), 
the  elapsed  time,  the  grid  size,  the  viewing  azimuth,  and  the  viewing  elevation. 


7.2.  Display  Output  Files 

Display  will  create  one  file  that  is  readable  by  Analyzer2  [6,7],  see,  and  see2.  This  file  contains  a 
temperature  for  each  pixel  of  the  file  just  ray-traced.  The  first  line  of  this  file  contains  two  integers  that  are 
the  width  and  height  of  the  grid.  They  should  be  the  same  since  display  uses  a  square  grid.  Each 
succeeding  line  will  contain  a  floating  point  number  that  is  a  temperature  of  a  pixel.  The  first  temperature 
is  for  the  top,  left  pixel.  The  pixels  move  to  the  right  and  then  down. 


7,3.  Display  Example 

This  example  again  uses  the  mged  file  test.veh.g.  Note  that  only  one  object,  vehicle,  was  used.  The  air 
regions  were  removed  from  the  PRISM  file;  therefore,  the  air  should  not  be  ray-traced. 
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$  display  test.veh.g  vehicle 

Type  of  output  file  to  be  read  0=  >PRISM,  1=  >  generic. 

0 

Enter  name  of  the  PRISM  output  file  to  be  read  (26  char  max). 
veh.prm 

Enter  the  number  of  regions  in  the  PRISM  file,  must  be  more 
than  eight  (not  including  the  background). 

18 

Enter  name  of  region  #  &  name  file  to  be  read  (26  char  max). 
veh.f.rnn 

Enter  name  of  output  file  (26  char  max). 
veh.dis 

Enter  the  elapsed  time  to  create  graphical  representation  of 
4 

Number  of  regions  (including  the  background):  19 
7/19/1984  7.000000:0.000000 
Prism  owl  file  read. 

Region  #  &  name  file  opened. 

The  number  of  regions  read  from  the  output  file  and  the  region  #  &  name 
file  was  the  same,  18  (does  not  include  background  in  number). 

Building  directory. 

File:  test.veh.g 

Database  Title:  test  vehicle  for  use  with  irprep  programs 
vehicle  loaded. 

The  number  of  regions  read  from  the  output 
file,  the  region  #  &  name  file,  and  the  .g 
file  are  all  equal.  The  number  of  regions 
read,  including  the  background  is  19 
Preparation  started. 

Minimum  &  maximum  X:  -8001.000000  -  8714-000000 

Minimum  &  maximum  Y:  -1850.000000  -  1850.000000 

Minimum  &  maximum  Z:  0.000000  -  2500.000000 

Center  of  bounding  sphere:  856.500000,  0.000000,  1250.000000 

Radius  of  bounding  sphere:  8829.551456 

Enter  multiplication  factor  for  radius. 

.75 

Enter  grid  size. 

512 

Enter  azimuth  &  elevation. 

35  25 

gridsize:  512  x  512 
azimuth:  85.000000  degrees 
elevation:  25.000000  degrees 
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The  command  line  for  using  display  is 
display  model. g  objects... 

where  model. g  is  the  mged  geometry  (test.veh.g)  and  objects  are  all  the  groups  except  air  (vehicle).  A  file 
needs  to  be  read  that  contains  the  temperature  information.  This  may  be  a  PRISM  or  generic  file.  In  this 
case  a  PRISM  file  (veh.prm)  is  read.  When  a  PRISM  output  temperature  file  is  used  the  number  of  regions 
needs  to  be  entered  since  this  is  not  contained  in  the  PRISM  temperature  output  file.  The  user  will  also 
need  to  enter  the  region  number  and  name  file  that  was  created  by  firpass.  The  user  than  needs  to  enter 
the  name  of  the  file  to  be  created  (veh.dis)  and  what  the  elapsed  time  is  for  ray-tracing.  Once  this  has  been 
done  information  about  the  file  is  printed  out.  The  user  is  then  asked  to  enter  a  multiplication  factor.  This 
makes  the  vehicle  smaller  or  larger.  Often  1.0  will  work,  but  this  is  something  the  user  needs  to  experiment 
with.  .Next  the  user  needs  to  enter  the  grid  size  for  ray-tracing,  256  and  512  are  good  choices.  At  this  time 
see  and  8ee2  which  actually  display  the  image  are  not  set  up  to  handle  any  file  larger  than  512.  The  user 
al.so  needs  to  enter  the  viewing  angles,  azimuth  and  elevation.  Now  display  will  create  a  file  readable  by 
Analyzer2,  see,  and  8ee2.  This  file  does  not  appear  in  an  appendix  due  to  its  length.  If  a  grid  spacing  of 
512  is  used  there  are  262145  lines  of  data. 


8.  SEE 

8.1 .  See  Input  File 

See  takes  advantage  of  X-Windows;  therefore,  when  running  see  the  user  must  be  using  X-Windows.  See 
reads  the  file  that  was  created  by  display  and  draws  a  picture  on  the  screen  giving  an  appropriate  color  to 
eacli  pixel  based  on  the  temperature  of  that  pixel.  In  order  to  run  see  the  user  must  have  a  compiled  copy 
of  see  and  a  file  created  by  display. 


8.2.  See  Output 

See  generates  an  image  in  one  of  the  four  color  scales  listed  below. 

gray 

red 

black-blue-cyan-green-yellow-white 

black-blue-magenta-red-yellow-white 

8.3.  See  Example 
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See  is  run  by  typing  see  (while  in  an  X-Windows  environment)  and  answering  the  questions  as  they  are 
asked. 

$  see 

Enter  name  of  file  to  be  read  (26  char  max). 
veh.dis 

Indicate  type  of  color  shading  to  use. 

0  -  gray 

1  -  red 

2  -  black-blue-cyan-green-yellow-white 
S  -  black-blue-magenta-red-yellow-white 

0 

Do  you  wish  to  create  a  piz  file  (0-no,  1-yes)? 

0 

Zeroing  color  info  array  -  finished  zeroing 
Setting  up  color  scale  -shades  of  gray  -  finished. 

Reading  file  -  file  read. 

Width:  512 
Height:  512 
Finding  min  &  max. 

Minimum:  22.590000 
Maximum:  ^6.830000 
Finding  pixel  bins  •  found  pixel  bins. 

Putting  color  info  in  arrays  -  color  info  in  arrays. 

$ 

First  the  user  enters  the  name  of  the  file  created  by  display  (veh.dis)  and  then  which  color  scale  is  to  be 
used.  The  user  is  asked  if  a  pix  file  should  be  created.  If  1  (yes)  is  entered  the  user  is  prompted  for  the  name 
of  the  pix  file.  This  file  can  be  used  in  conjunction  with  other  routines  in  the  BRL-CAD  package.  Next  the 
program  prints  out  statements  indicating  the  stage  of  the  program.  It  will  take  a  little  time  to  read  the  file. 
Then  the  image  is  displayed  on  the  screen  (Figure  2). 


9.  SEE2 

9.1.  See2  Input  File 

The  difference  between  see  and  8ee2  is  that  8ee2  takes  advantage  of  the  SGI  graphics  instead  of  X- 
Windows.  Again,  this  program  uses  the  file  created  by  display  and  creates  an  image  using  one  of  the  color 
scales  listed  below. 


gray 

black-blue-cyan-green-yellow-white 
black-blue- magenta-red-yellow- white 
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The  user  will  need  a  compiled  copy  of  see2  and  the  file  created  by  display. 


9.2.  Se«2  Output 

The  output  of  8ee2  is  an  image  in  one  of  the  color  scales  listed  above. 

9.3.  See2  Example 

See2  is  run  by  typing  see2  (on  an  SGI)  and  answering  the  questions. 

$  see2 

Enter  name  of  file  to  be  rend  (26  char  max) 
veh.dis 

Indicate  color  scale  to  be  used. 

0  -  gray 

1  -  black-blue-cyan-green-yellow-white. 

2  -  black-blue-magenta-red-yellow-white. 

0 

Print  scan  line  number  (0-yes,  1-no)? 

1 

Do  you  wish  to  create  a  pix  file  (0-no,  1-yes)? 

0 

Setting  color  scale  -  gray  scale  -  set. 

Reading  file  -  file  read. 

Width:  512 
Height:  512 
Finding  min  &  max. 

Minimum:  22.590000 
Maximum:  46.8S0000 
Finding  pixel  bins  -  pixel  bins  found. 

Setting  color  for  each  pixel  -  colors  found. 

Press  'x'  return  to  end.  z 

THE  END 

$ 

First  the  user  enters  the  name  of  the  file  created  by  display  (veh.dis)  and  then  selects  a  color  scale  (gray). 
Then  the  user  is  asked  if  the  scan  lines  should  be  printed,  the  answer  should  probably  be  no  since  printing 
scan  lines  slows  the  proce.ss  down.  The  printing  of  the  scan  lines  was  only  added  for  diagnostic  purposes. 
The  user  then  entered  0  to  indicate  a  pix  file  was  not  to  be  created.  If  the  user  had  entered  1  there  would 
have  been  a  prompt  to  enter  the  pix  file  name.  The  image  will  be  displayed  and  the  user  presses  z  and 
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return  to  end  the  process.  See  figure  3  for  the  image  created  by  8ee2. 


10.  PICTX  AND  PICTSGI 

Pictx  and  pictsgi  combine  display  and  sec  or  8ee2  together  into  one  program.  Pictx  combines  display 
and  see  and  needs  to  be  run  in  an  X-Windows  environment.  Pictsgi  combines  display  and  see2  and  runs 
on  an  SGI.  Pictx  simply  calls  display  and  see  while  pictsgi  calls  display  and  see2.  Start  pictx  by 
typing  pictx  and  pictsgi  by  typing  pictsgi.  The  user  only  needs  to  answer  the  questions  as  they  appear. 
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E]  Coior  Scale 


Figure  3.  Image  Displayed  by  See2 
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APPENDIX  A  -  Materials  Properties  File 
(See  section  3.1  for  a  description  of  this  file.) 
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absorpt iv ity 
emissivlty 

thermal  conductivity  (W/raK) 
material 


APPENDIX  B  -  Output  Files  from  firpass 
(See  section  3.3  for  a  description  of  these  files.) 
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2  Facet  file  for  use  with  PRISM. 

FN  DESCRIPTION  TY  AREA  MASS  SPHEAT  El  E2  ABSOR 

SN(X)  SN(Y)  SN(Z)  CONY  K  L  SHAPE 

1  /vehicle/hull/r.hullOl  1  6.320  939.456  460.200  0.940  0.940  0.600 

0.000  1.000  0.000  1.000  0  0  0.000  1.176  0.000 
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APPENDIX  C  -  Output  Files  from  secpass 
(See  section  4.3  for  a  description  of  these  files.) 
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